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; I.  INTRODUCTION 

It  seems  self-evident  that  the  use  of  aciiial  measurements  of 

( aerodynamic  drag  on  a satellite  rather  than  a force  model  based  on  an 

t 

' ■ atmospheric  model  would  provide  improved  capability  for  post-flight 

I 

■ orbit  reconstruction  and  near-term  ephemeris  prediction.  However,  in 

I order  to  justify  the  use  of  drag-measuring  devices  on  an  operational  basis, 

[ 

f a decision  would  have  to  be  made  as  to  whether  the  improved  capability 

i 

was  sufficiently  great  to  warrant  the  additional  expense.  This  means,  of 

[ 

I course,  that  one  must  be  able  to  define  what  constitutes  a ''sufficient'' 

improvement  and  also  have  quantitative  data  upon  which  to  base  a decision. 

The  purpose  of  the  study  reported  here  was  to  obtain  such  quanti- 
tative information.  Realistic  conditions  were  achieved  by  using  calori- 
meter data  from  an  actual  flight.  However,  in  order  not  to  obscure  the 
. interpretation  of  the  results,  it  was  necessary  to  assume  that  the  calori- 

meter data  were  "perfect,  " i.  e.  that  any  errors  in  orbit  fitting  or  pre- 
dicting were  due  to  factors  other  than  drag-measurement  errors.  In  this 

' sense,  the  results  of  the  study  are  indicative  of  the  "best"  possible  per- 

formance that  should  be  expected  and  that  in  ac‘:ual  practice  one  would 
have  to  allow  for  errors  of  measurement. 

This  study  proceeded  along  two  -^er.eral  lines:  (1)  simulated  orbit 
calculations  using  the  TRACE  program,  and  (2)  evaluation  of  the  relation- 
ship between  prediction  errors  and  the  drag  model.  The  latter  effort  in- 
cluded an  analysis  of  the  drag  data  in  search  of  systematic  patterns  of 
behavior  that  might  be  used  to  improve  predictive  capability. 

-9- 
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In  Section  II  we  shall  describe  in  detail  the  objectives  and  logistics 
of  the  orbital  calculations,  while  in  Section  III  we  will  present  the  results 
and  discuss  the  relationship  between  orbit  prediction  errors  and  the  pre- 


j dictive  drag  model.  A summary  of  the  findings  along  with  reservations  |1 

I and  recommendations  will  be  presented  in  Section  IV.  ! 

t I 

f 

I I 

! i 


II.  ORBIT  SIMULATION  STUDIES 

1 . Objective  and  Method 

The  basic  objective  of  this  phase  of  the  study  was  to  determine  the 
errors  in  3-rev  predictions  that  would  typically  result  from  the  use  of  in 
situ  drag  measurements  as  opposed  to  the  operational  approach  which 
relies  on  the  'DENSER  model  to  define  the  drag  force.  It  is  assumed  that 
the  drag  data  (in  this  case  calorimeter  measurements)  are  available  only 
in  the  6-rev  period  immediately  preceding  the  3-rev  predict  span  and  that 
these  data  may  be  used  to  derive  (1)  the  vehicle  position  and  velocity 
vector  at  the  start  of  the  predict  span,  and  (2)  a drag-force  model  for 
the  predict  span. 

To  obtain  the  "error"  one  must,  of  course,  have  a "reference" 
trajectory  representing  the  "true"  world.  The  reference  trajectory  was 
specified  by  using  9-revs  of  calorimeter  data  to  represent  the  drag  force 
along  with  a nominal  geopotential  model  and  initial  conditions  suitable  for 
a low -altitude  orbit. 

The  sequence  of  TRACE  calculations  in  shown  in  Fig.  1.  First, 
a 9-rev  reference  trajectory  was  generated  by  TRACE.  Artificial  tracking 
data  with  random  noise  added  were  generated  during  the  first  6-revs.  The 
tracking  station  locations,  data  type  and  rate,  sigmas,  etc.  were  taken 
to  be  representative  of  the  ACES  network. 

The  second  step  in  the  procedure  was  to  use  the  simulated  tracking 
data  to  obtain  a 7-parameter  orbit  fit  over  the  first  6-revs.  In  the  flow- 
chart, tj  and  t^  are  the  start  and  stop  times  of  the  fit  span.  Orbit  fit- 
ting was  done  using  'DENSER  (the  operational  model)  and  also  using 


» 

i 

i‘ 
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calorimeter  data.  The  fitted  orbits  were  differenced  with  respect  to  the 
reference  orbit,  yielding  radial,  intrack,  and  crosstrack  errors.  (It  is 
to  be  noted  that  since  the  same  calorimeter  data  are  used  to  gene  ate  the 
reference  trajectory  and  obtain  an  orbit  fit,  the  resulting  fit  necessarily 
represents  the  best  possible  fit  one  might  realistically  hope  for.  The 
only  factor  that  precludes  a ''perfect''  fit  is  the  random  noise  in  the  track- 
ing data.  ) 

The  third  step  in  the  procedure  involved  predicting  ahead  for  3-revs 
(time  t^  to  t^  in  the  flowchart  of  Fig.  1 ) using  the  fitted  vector  at  t^ 
to  initialize  the  predictions.  The  predicted  and  reference  trajectories 
were  differenced.  Several  different  drag  models  were  tested  in  the  predict 
span:  (1)  'DENSEL,  (2)  a model  that  represented  the  average  of  all  6-revs 
of  calorimeter  data  in  the  fit  span,  called  Model  B,  and  (3)  a model 
that  represented  the  calorimeter  data  of  only  the  last  rev  of  the  fit  span, 
called  Model  C.  Models  B and  C were  actually  smoothed  and  edited 
versions  of  the  raw  calorimeter  data.  The  manner  by  which  the  calori- 
meter data  were  smoothed  and  then  transformed  to  working  density  models 
is  described  in  detail  in  a preliminary  report  (Ref.  1 ). 

In  the  calorimeter  data  base  that  was  available  for  our  use,  only 
5 sets  of  9 consecutive  orbits  of  data  were  found.  The  study  reported 
here  is  based  on  only  those  5 sets.  In  an  operational  sense  it  would  be 
quite  simple  to  devise  a routine  to  handle  cases  of  missing  data. 

For  each  of  the  5 sets  of  data  the  following  combinations  of  drag 
models  were  used  to  fit  and  predict: 


INPUT 

DRAG  MODEL 


CALORIMETER 

DATA 


REFERENCE  TRAJECTORY 
time  t,  to  t» 


ORBIT 

DIFFERENCING 


Simulated 
Tracking  data 


OUTPUT 


' DENSE  L, 
CALORIMETER 


ORBIT  FITTING 


Initial 
Conditiona 
At  t. 


ORBIT 

DIFFERENCING 


OUTPUT 


'DENSEL, 

:AL0RIMETER 


ORBIT  prediction! 


1.  Flow  Chart  of  Orbit  Simulation  Tests 


1. 


'DENSEL;  'DENSEL  ( operational  method  ) 

2.  Calorimeter  data;  Model  B 

3.  Calorimeter  data;  Model  C 

4.  Calorimeter  data;  'DENSEL. 

In  a few  instances  other  combinations  of  models  were  tested.  Such  cases 
will  be  identified  and  discussed  in  the  next  section.  Whenever  the  'DENSEL 
model  was  used  in  prediction,  the  ballistic  factor  ( C^^A/W  ) was  scaled 
according  to  the  drag  parameter  obtained  in  the  orbit  fit  with  'DENSEL. 

2.  Orbit  Characteristics 

Initial  conditions  for  the  reference  trajectories  were  selected  so 
as  to  reconstruct  as  closely  as  possible  the  actual  orbits  corresponding  to 
the  calorimeter  data.  Briefly,  the  orbits  were  near-circular  and  low 
altitude  (roughly  68  by  200  nmi  ),  sun- synchronous  ( i ~ 96°  ),  and  with 
perigee  at  middle  latitudes.  In  one  set  of  calculations  we  tested  two 
different  initial  values  of  the  right  ascension  of  the  ascending  .lode.  This 
had  the  effect  of  calling  different  tracking  stations  into  play.  The  results 
of  this  change  were  insignificant  relative  to  other  effects  being  studied 
and,  therefore,  will  not  be  given  further  consideration  in  this  report. 

The  orbit  numbers  of  the  data  samples  and  other  pertinent  infor- 
mation are  given  in  Table  1. 
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TABLE  1.  Orbit  Characteristics  of  the  Data  Samples 


Sample 

Number 

Orbit 

Numbers 

Perigee 

Height  (nm) 

Perigee 
Latitude  (deg) 

Apogee 
Height  (nm) 

1 

179-187 

69 

42 

217 

2 

325-333 

69 

46 

210 

3 

502-510 

69 

44 

208 

4 

761-769 

72 

52 

205 

5 

792-800 

72 

56 

206 

3. 


Calorimeter  Data  Base 


The  data  base  for  this  study  consisted  of  calorimeter  measurements 
of  the  aerodynamic  heating  on  a low-altitude  satellite  obtained  over  a period 
of  slightly  less  than  2 months.  Roughly  300  data  passes  were  acquired  during 
this  interval,  but  only  5 sets  of  9 consecutive  orbits  contained  sufficient  cover- 
age  about  perigee  to  be  useful  for  the  orbit -simulation  calculations.  However, 
all  useful  passes  of  data  were  included  in  studies  related  to  atmospheric  be- 
hehavior.  An  example  of  a fitted  pass  of  data  is  shown  in  Fig.  2. 

The  fundamental  premise  of  this  study  is  that  the  atmospheric  density 
derived  from  the  heat  flux  data  is  a "perfect"  representation  of  the  true  atmo- 
sphere. This  assumption  is,  of  course,  not  very  realistic  since  our  know- 
ledge of  the  vehicle  drag  coefficient  and  the  surface  thermal  accommodation 
coefficient  (along  with  their  altitude  dependences)  is  slightly  less  than  perfect, 
l^^rther,  uncertainties  associated  with  the  nneasurements  were  unknown  to  us. 
These  factors,  therefore,  must  be  considered  to  be  outside  the  scope  of  the 
present  study,  and  the  results  to  be  presented  necessarily  reflect  optimum 
conditions. 
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Fig.  2.  Example  of  Calorimeter 
Data  Fitting 
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III.  RESULTS  AND  INTERPRETATION 

1 . TRACE  Calculations 

Graphs  of  the  intrack  and  radial  residuals  for  the  five  data  sets  are 
presented  in  Appendix  A.  For  each  data  set  there  are  four  combinations 
of  fit-and-predict  drag  models  as  discussed  previously.  In  the  graphs  the 
thick  vertical  line  separates  the  fit  and  predict  intervals;  arrows  along  the 
lower  horizontal  axis  indicate  the  perigee  positions  within  the  predict  span. 

Root-mean-square  errors  in  the  predict  span  are  presented  in 
Table  2.  Errors  at  perigee  are  given  in  Tables  3a  and  3b  for  the  radial 
and  intrack  components,  respectively.  Without  exception,  orbit-fitting 
with  calorimeter  data  resulted  in  better  fits  than  the  operational  model 
and,  thus,  provided  better  starting  vectors  for  the  predict  spans,  RMS- 
errors  of  the  radial  component  in  the  predict  span  were  typically  a factor 
of  2 or  better  when  calorimeter  data  were  used;  the  improvement  in  the 
RMS-error  of  the  intrack  component  was  typically  a factor  of  3 or  more. 

Focusing  on  the  errors  at  perigee  (Tables  3a,  3b)  we  can  make  the 
following  observations; 

1.  With  but  one  exception  (fit  with  calorimeter  data,  predict  with 
MODC),  the  average  absolute  radial  error  at  perigee  was  less 
when  calorimeter  data  were  used  in  the  fit  than  when  'DENSE L 
model  was  used. 

2.  The  average  absolute  radial  error  at  perigee  was  significantly 
less  than  the  RMS-error  over  the  predict  span  regardless  of 
what  model  was  used  to  fit  or  predict. 
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TABLE  Z.  RMS  Errors  (in  feet)  in  the  Predict  Span 


Predict  Span 
Rev.  No. 


185  - 187 

radial 

intrack 


331  - 333 

radial 

intrack 


508  - 510 

radial 

intrack 


767  - 769 

radial 

intrack 


798  - 800 

radial 

intrack 


*CALCR/  CALOR/ 
MOD  B MOD  C 


288 

2828 


CALOR/ 

'DENSEL 


'DENSEL/ 

•DENSEL 


10 

258 

n 

3739 

434 

4848 


1107 
14, 720 


* CALOR /MODE  means  calorimeter  data  were  used  in  the  fit 
span.  Model  B in  the  predict  span. 


Predict  Span 
Rev.  No. 

CAL OR/ 
MOD  B 

CAL OR/ 
MOD  C 

CALOR/ 

'DENSEL 

'DENSEL/ 

'DENSEL 

185 

-14 

1 

0 

186 

1906 

-166 

1092 

187 

4942 

-1331 

24V9 

9153 

Abs.  Ave. 

2287 

499 

1190 

5877 

331 

-29 

-28 

-16 

1446 

332 

-326 

362 

768 

3585 

333 

-1480 

576 

1792 

5971 

Abs.  Ave. 

612 

322 

859 

3667 

508 

-11 

-4 

-10 

4056 

509 

2565 

1430 

3684 

12, 125 

510 

4430 

1022 

7798 

22, 384 

Abs.  Ave. 

2335 

812 

3831 

12, 855 

767 

-27 

-24 

-25 

499 

768 

810 

-316 

367 

1656 

769 

2354 

-1066 

1056 

3126 

Abs.  Ave. 

1063 

469 

483 

1760 

798 

-5 

7 

-5 

678 

799 

-224 

-156 

69 

1628 

800 

-951 

-738 

-65 

2741 

Abs.  Ave. 

393 

300 

46 

1682 
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The  absolute  radial  error  at  perigee  did  not  consistently  increase 
with  time  in  the  predict  span. 


i 


I 


4.  The  average  absolute  intrack  error  at  perigee  was  always  less 
using  calorimeter  data  to  fit  than  using  the  'DENSEL  model. 

5.  The  average  absolute  intrack  error  at  perigee  did  not  differ 
significantly  from  the  RMS-error  over  the  predict  span. 


I 


6.  The  absolute  intrack  error  at  perigee  increased  with  time  in 
the  predict  span. 

7.  Radial  errors  at  perigee  did  not  correlate  with  any  of  the  3 drag 
models  used  in  prediction.  By  constrast,  the  intrack  errors  at 
perigee  were  always  the  least  for  Model  C (based  on  the  last  rev 
of  calorimeter  data  in  the  fit  span). 

To  provide  better  insight  on  the  relationship  between  errors  and  the 
drag  prediction  model,  we  have  plotted  the  RMS-errors  and  the  perigee  errors 
of  both  the  radial  and  intrack  components  against  the  following  quantities: 

(1)  the  initial  error  of  the  predict  span,  (2)  the  mean  level  of  tl  - true 
density  in  the  predict  span,  and  (3)  the  density  variability  withing  the  predict 
span.  In  order  to  avoid  ambiguities  in  interpretation,  the  drag  model  used 
for  prediction  in  all  cases  was  the  operational  model  'DENSEL. 


Figure  3a  shows  the  dependence  of  the  RMS-radial  error  in  the  predict 
span  on  the  initial  radial  error  in  the  starting  vector.  The  5 circles  on  the 
graph  represent  the  5 different  data  sets  that  were  fit  with  calorimeter  data; 


J. 
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the  X's  represent  the  fits  with  'DENSEL.  As  mentioned  previously,  all 
predictions  were  made  with  'DENSEL  and  the  ballistic  coefficient  was 
scaled  according  to  the  fitted  drag  parameter.  On  the  scale  used  for  the 


abscissa,  the  initial  errors  in  the  radial  component  | ARq  | for  the 
calorimeter  fits  are  indistinguishable  and  appear  to  be  near-zero.  The 
actual  values  ranged  from  a few  tenths  of  a foot  to  a few  feet.  Since, 
however,  the  RMS-errors  were  uncorrelated  with  | AR^  [ for  the  calori- 
meter cases,  we  saw  no  advantage  in  expanding  the  scale  near  zero.  In 
the  cases  where  'DENSEL  was  used  to  fit  the  tracking  data,  there  is  a well- 
defined  relationship  between  the  initial  error  and  the  RMS-error.  If  we 
were  to  average  calorimeter  points,  that  average  ( ~ 200  ft  ) would  fit  in 
nicely  with  the  other  points.  Essentially  the  same  remarks  may  be  applied 
to  the  intrack  component  shown  in  Fig.  3b. 

In  Figs.  4a  and  4b  the  average  of  the  absolute  errors  at  perigee  in 
the  predict  span  are  plotted  against  the  initial  errors  in  the  radial  and  intrack 
components,  respectively.  Again,  the  initial  errors  corresponding  to  the 
calorimeter  data  fits  are  so  much  less  than  the  errors  from  the  'DENSEL 
fits  that  their  points  are  crowded  near  zero.  As  before,  expansion  of  the 
scale  near  zero  would  serve  no  purpose  since  the  errors  appeared  to  be  un- 
correlated. As  opposed  to  Figs.  3a  and  3b,  there  is  no  well-behaved  re- 
lationship between  the  perigee  errors  and  the  initial  errors  in  the  radial 
component.  The  intrack  component  does  exhibit  a dependence,  but  with  some 
scatter. 
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In  Figs.  5a  and  5b  we  show  the  RMS-crrors  in  the  predict  interval 


plotted  against  the  mean  density  of  the  predict  interval.  The  calorimeter 
and  the  'DENSEL  points  both  show  that  the  RVS-error  increases  with  the 
mean  le\  el  of  the  density,  which  is  intuiti\  ely  what  one  would  expect. 

With  regard  to  the  perigee  errors,  it  is  seen  in  Figs.  6a  and  6b 
that  the  intrack  component  does  increase  with  the  density  while  the  radial 
component  does  not. 

In  Figs.  7a.  b and  8a,  b we  show  the  RMS  and  perigee  errors 
against  the  density  \ariability  within  the  predict  interval,  the  latter  being 
defined  as  the  RMS -de\iation  from  the  mean.  In  all  cases  we  obserxe  a 
definite  tendency  for  errors  to  grow  with  density  variability. 

To  summarize  Figs.  3-8,  we  note  that  intuiti\ely  one  would  expect 
errors  to  grow  in  all  cases,  i.e.  , with  initial  error,  with  mean  density, 
and  with  \ariability  in  the  dt'nsity.  Indeed,  this  was  \'erified  for  the  RMS- 
error,  hich  e s sentially  is  an  average  o\er  the  entire  predict  span,  and  the 
intrack  error  at  perigee,  but  not  the  radial  error  at  perigee.  In  the  latter 
case  only  the  density  variability  factor  appeared  to  influence  the  radial 
error  at  perigee  in  a systematic  way. 

These  orbit  simulation  studies  indicate  that  both  the  starting  \ector 
of  the  predict  span  and  the  density  within  the  predict  span  are  impe^rtant  factors 
for  goc:>d  ephemeris  prediction.  Furthc-r,  it  is  not  only  the  mean  density 
that  is  important,  but  time -dependent  clianges  as  well. 


\ 
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/.tmospheric  Studies 


At  the  low  perigee  altitudes  of  the  orbits  under  study,  very 

little  is  known  of  the  density  variations  that  occur  on  a time  scale  of  a 

few  hours  except  that  they  do  exist  and  can  be  substantial  ( - 10%  of  the 

mean  level).  At  higher  altitudes  ( 100  nm  ) it  is  known  that  such  short- 

period  \ariations  tend  to  correlate  with  magnetic  activity.  With  the  pres- 

, ent  set  of  calorimeter  data  we  found  the  short-term  orbit-inte" rated 

! 

density  variations  to  be  rather  poorly  correlated  with  Kp  fluctuations, 
although  the  level  of  magnetic  activity  throughout  the  data  interval  w'as 
relatively  low.  How'ever,  on  a longer  time  scale  of  a day,  the  orbit- 
integrated  density  did  show  a relationship  with  daily  magnetic  activity 
as  illustrated  in  Fig.  9.  (In  this  and  following  figures  the  density  has 
been  adjusted  to  compensate  for  changes  related  to  altitude.  This  nor- 
malization procedure,  w'hich  is  frequently  used  in  aeronomical  studies, 
is  discussed  in  Appendix  B.  ) 

Since  orbit-to-orbit  density  changes  could  be  associated  with 
spatial  as  well  as  temporal  variations,  we  investigated  the  possibiiiiy 
of  geographically-related  features  in  the  density  distribution.  A first- 
cut  analysis  involved  plotting  the  orbit-integrated  density  vs  longitude 
of  minimum  altitude;  this  is  shown  in  Fig.  10.  Three  horizontal  lines  are 
shown  for  referen  e:  the  upper  line  delineates  the  upper  quartile  of  the 
data,  the  middle  line  is  the  mean  (also  the  median),  and  the  lower  line 
delineates  the  lower  quartile.  The  machine-plot  is  in  a sense  misleading 
since,  if  more  than  one  point  falls  within  a specified  grid,  only  a single 
asterisk  is  shown.  Thus,  the  number  of  points  that  should  lie  between 
the  upper  and  lower  lines  is  greater  than  is  actually  shown. 
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Ly  Average  o£  Orbit-lntegralecl  Density  (eorreelecl  fer  rlille  renees 
erigee  lieigli!)  and  Daily  Snni  o£  Kp 


The  plot  indicates  that  the  density  near  72°E  longitude  is  greater 
by  about  10%  than  the  density  near  20°E.  There  appears  to  be  a minimum 
near  20°E  longitude,  but  this  is  speculative  since  there  are  too  few  points 
to  the  west  of  the  apparent  minimum.  The  density  in  the  sector  from  144°E 
to  about  240°E  also  appears  to  be  lower  than  average;  however,  again  there 
are  too  few  points  to  be  certain  of  this  feature. 

Plotting  the  density  against  both  latitude  and  longitude  of  minimum 
altitude  as  in  Fig.  11  fails  to  elucidate  matters.  Here,  circles  represent 
the  upper-most  quartile;  asterisks  are  used  for  the  upper-middle  quartile, 
plus -signs  for  the  lower -middle  quartile,  and  minus- signs  for  the  lower 
quartile.  From  about  35  N to  50°N,  where  the  majority  of  points  exist, 
the  symbols  are  more  or  less  distributed  in  buck-shot  fashion.  There  is 
a definite  lack  of  upper -quartile  symbols  in  the  144-240°E  sector,  but 
as  mentioned  above  the  data  coverage  there  was  very  sketchy.  Similarly, 
there  appears  to  be  a zone  of  high  density  between  20°  and  24°N  and  a 
zone  of  low  density  between  24°  and  30°N:  but  again  the  poor  data  coverage 
may  be  biasing  the  statistics.  In  all,  therefore,  where  there  was  good  data 
coverage,  there  was  no  well-defined  geographically-related  structure. 
Possible  regions  jf  "highs"  or  "lows"  were  compromised  by  the  paucity 
of  data. 

Although  no  outstanding  geographical  or  magnetic -index  related 
effects  were  found,  it  was  determined  from  the  calorimeter  data  base  that 
the  density  at  the  heights  in  question  did  exhibit  a certain  degree  of  what 
might  be  called  persistence.  That  is,  given  the  orbit-integrated  density  for 
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Distribution  of  Altitude -Corrected  Orbit-Integrated  Density  in  Latitude 
and  Longitude  of  Perigee:  0 = upper  quartile;  -i'  = upper-middle 
quartile;  + = lower-middle  quartile;  - = lower  quartile. 


a particular  orbit,  we  found  that  in  68  of  155  cases  ( 44%  ) the  integrated 
density  of  the  following  orbit  did  not  differ  by  more  than  + 5%.  Further,  in 
88%  of  the  cases,  the  integrated  density  was  within  + 15%  of  that  of  the 
preceding  orbit.  Similarly,  in  ~ 30%  of  the  cases  available  for  study,  the 
density  of  the  second  and  third  orbits  following  the  given  orbit  did  not 
change  by  more  than  ^ 5%.  This  "persistence"  feature  explains  why  in 
nearly  all  our  calculations  the  prediction  Model  C yielded  the  lowest  errors. 
Model  C,  it  is  recalled,  was  simply  the  measured  density  of  the  last  orbit 
in  the  fit  span. 

3.  Further  Considerations  of  Short-Term  Prediction 

The  question  of  the  relative  importance  of  the  starting  vector  as 
opposed  to  the  accuracy  of  the  drag  model  in  the  predict  span  was  investigated 
by  comparing  several  TRACE  runs  having  different  combinations  of  fit/predict 
drag  models.  The  comparison  was  carried  out  for  only  one  of  the  five  data 
samples,  and  so  the  results  must  be  viewed  and  interpreted  with  a fair  amount 
of  discretion. 

Basically,  the  purpose  was  to  compare  "perfect"  with  "operational" 
drag  force  models,  where  the  "perfect"  model  was  the  calorimeter  data  that 
were  used  to  generate  the  reference  trajectory,  and  the  operational  model  was, 
of  course,  'DENSEL.  The  results  are  presented  in  Table  4.  Both  RMS- 
errors  and  perigee  errors  in  the  predict  span  are  tabulated.  The  data  sample 
used  was  Revs  761-769. 

For  this  particular  sample,  having  either  a perfect  starting  vector  or 
a perfect  drag  model  yielded  essentially  similar  improvements  in  RMS-errors 
as  compared  with  the  operational  approach.  On  the  other  hand,  the  accuracy 
of  the  starting  vector  appears  to  be  important  insofar  as  perigee  errors  are 


concerned. 
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As  mentioned  before,  we  must  be  careful  not  to  over -interpret  the 
results  because  we  are  dealing  with  a single  data  set.  A cautious  conclusion 
that  may  be  drawn  from  this  exercise  is  that  both  the  starting  vector  and  the 
drag  model  are  important,  and  thus  future  endeavors  to  improve  short-term 
orbit  prediction  should  not  overlook  either  of  these  factors. 


Table  4.  RMS-Errors  and  Mean  Absolute  Errors  at  Perigee 
(in  feet)  in  the  Predict  Span  of  Data  Set  4 ( Revs  767-769  ) 


'DENSEL/'DENSEL 

CALCR/'DENSEL 

' DENSER /CALOR 

Radial 

RMS 

470 

250 

257 

Perigee 

147 

58 

173 

Intrack 

RMS 

3086 

1194 

1960 

Perigee 

1760 

483 

1286 

An  interesting  aspect  of  prediction  errors  that  can  be  seen  in  the  figures 
in  Appendix  A is  the  complexity  and  sensitivity  of  the  time -dependent  behavior 
related  to  different  predictive  drag  models.  That  is,  given  a specific  starting 
vector,  the  drag  model  affects  not  only  the  envelope  of  the  error  oscillations, 
but  the  phase  as  well.  The  source  of  this  complex  response  can  be  understood 
at  least  qualitatively  by  examining  simplified  analytical  expressions  for  the 
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orbital  perturbations  due  to  air  drag.  Specifically,  the  radial  distance  r 
at  a position  on  the  orbit  where  the  true  anomaly  is  0 is  related  to  the 
distance  r^  for  a drag-free  orbit  by  (Ref.  Z,  eqn. 


1^ 
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1_ 

r 
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/' 


1 - cos  ( ^ - 0 ) 


2 1/2 

1 ( 1 + Ze  cos  e + e ) ' 

! 2. 

( 1 + e cos  0 ) 


P ( G ) de 


(1) 


where  p is  the  density,  e is  the  eccentricity,  and  6 is  a parameter 
similar  to  the  so-called  ballistic  parameter.  From  (1),  it  is  straight- 
forward to  derive  the  difference  in  radial  distance  after  one  revolution  o’  dng 
to  different  density  models: 


where 


G { e,  0 ) = 


0 + Ztt 
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- cos  ( 0 - 0 ) 
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■7  n ^ 2 1 /Z 

-r  Ze  cos  0 -t  e ) 

( 1 + e cos  0 


(2) 


(3) 


A sketch  of  the  behavior  of  the  function  G is  shown  in  Fig.  IZ.  The  important 
feature  is  that  the  phase  depends  on  0 , the  position  in  the  orbit  where  the 
radial  difference  is  to  be  evaluated.  In  essence,  the  difference  in  the  radial 
component  (eqn.  2)  is  the  orbit -integ rated  difference  in  density  weighted  by  the 
function  G. 

In  the  simple  case  where  the  ratio  of  the  densities  ( 1 constant 

along  the  orbit,  eqn.  (2)  may  be  written 
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In  deriving  (4)  we  assume  that  the  density  is  given  by 


P = Pp  exp  1 - ( r - r^  ) / H ] 


where  is  the  density  at  perigee,  r^  is  the  perigee  distance,  and  H the 

density  scale  height;  the  parameter  P in  (4)  is  r^e/H.  In  this  simple  case 
the  integral  is  a positive -definite  quantity,  and  thus  the  radial  difference  has 
the  same  sign  for  all  0 depending  on  the  ratio  It  is  apparent  from 

Fig.  IZ  and  eqn.  (4)  that  for  given  values  of  the  parameters  P^^Pi’  ^ 
and  p , the  difference  is  a minimum  at  perigee  (0=0)  and  a maximum  at 

pi 

apogee  (0  = it  ).  This,  of  course,  is  not  surprising  to  anyone  with  any  famili- 
arity with  orbital  analyses,  as  it  is  well  known  that  the  perigee  height  decays 
very  slowly  under  the  action  of  air  drag  relative  to  the  decay  of  apogee. 


In  the  general  case  represented  by  eqn.  (2),  the  radial  difference 
depends  on  the  variation  of  the  density  difference  around  the  orbit  weighted 
by  the  function  G.  It  should  be  noted  that  even  if  the  orbit-integrated  difference 
in  the  density  is  zero,  the  radial  difference  can  be  non-zero  because  of  the 
function  G,  It  is  obvious  that  under  general  conditions,  where  ( p^  ' P^  ^ 
varies  in  some  cor.  plex  fashion,  it  is  no  straightforward  matter  to  infer  the 
effect  on  tne  ra  bal  istance.  To  even  predict  the  sign  of  the  change  would 
require  k'ovvlei  e of  the  behavior  of  the  d nsity  difference  around  the  orbit. 
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SUMMARY  AND  RECOMMFIsDATIONS 


Subject  to  the  constraints  and  assuniptions  niontioncd  earlier,  the  major 
ults  of  this  study  are  as  follows: 

Given  the  tracking  coverage  of  the  AGES  network,  direct  r’'>ea_'.ircnients 
of  the  drag  on  a low -altitude  vehicle  allow  the  orbit  to  be  fit  with  far 
greater  accuracy  than  would  be  obtainable  using  .'he  operational  drag 
model.  RMS-errors  in  the  fit  span  of  the  radial, intrack,  and  cross- 
track components  were  ' 1 ft,  10  ft,  and  < 1 ft,  respectively, 

using  drag  data,  and  > 10  ft,  > 100  ft,  :■  10  ft  using  the  'DENSER  model. 

The  accurate  starting  vectors  obtained  with  the  drag  measurements 
significantly  improved  the  3-rev  prediction  capability.  Using  the 
'DENSER  model  to  predict,  we  achieved  an  average  reduction  in 
RMS  error  of  a factor  of  2 in  the  radial  component  and  a factor  of  4 
in  the  intrack  component.  The  mean  absolute  error  at  perigee  was 
reduced  by  factors  of  3 and  4,  respectively,  for  the  radial  and  intrack 
components. 

The  level  of  the  density  and  the  variability  of  the  density  within  the 
predict  span  affect  prediction  accuracy.  On  a short  term  basis 
(i.  e.  orbit  by  orbit),  we  found  no  outstanding  correlation  be  ' ’een  the 
density  and  the  often-used  Kp  index  of  magnetic  activity.  We  also 
failed  to  find  compelling  evidence  of  permanent  geographically-related 
features  in  the  density.  We  did  find,  however,  that  the  statistical 
"persistence"  of  the  lower  thermosphere  can  be  exploited  to  improve 
short  term  orbit  prediction. 
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As  emphasized  in  this  report,  the  results  must  be  interpreted  as 
"optimistic"  since  no  allowances  have  been  made  for  errors  in  drag  measure- 
ments, satellite  tracking  (other  than  random  errors),  or  other  external  force 
models  (i.  e.  the  geopotential).  The  next  obvious  test  would  be  to  use  real, 
not  simulated,  tracking  data  in  a comparative  study  of  drag  data  vs  opera- 
tional model. 

Although  the  emphasis  of  this  study  has  been  on  satellite -borne  mea- 
surements of  drag,  it  should  be  noted  that  improvement  in  the  starting  vector 
of  a prediction  -interval  might  also  be  achieved  by  improved  tracking  capability. 

A most  important  issue  is  the  need  to  improve  the  drag  model  predic- 
tive capability  on  the  short-term  time  scale  (i.  e.  an  orbital  period).  Although 
it  was  found  that  the  persistence  feature  of  the  atmosphere  could  be  exploited 
in  orbit  prediction,  success  is  to  be  expected  only  on  a statistical  basis.  It 
must  also  be  noted  that  the  present  data  base  was  obtained  during  a period 
devoid  of  major  magnetic  activity,  and  so  we  have  not  been  able  to  analyze  the 
effects  of  such  disturbances  on  orbit  prediction.  Large  changes  in  density 
taking  place  on  a short-time  scale  can  have  profound  effects  on  low-altitude 
orbits:  thus  it  is  important  to  be  able  to  predict  the  occurrence  and  magnitude 
of  such  events.  A break-through  in  this  area  does  not  appear  to  be  imminent, 
but  since  the  drag  model  may  become  the  ultimate  limiting  factor  in  orbit 
prediction,  the  task  should  not  be  readily  abandoned. 
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APPENDIX  A.  Plots  of  Orbit  Differences 
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Plots  of  the  radial  and  intrack  errors  are  shown  on  the  following  pages. 
Figures  of  the  crosstrack  error  have  not  been  included  since  this  component 
is  generally  not  a source  of  significant  error. 

The  drag  models  used  in  the  fit  and  predict  intervals  and  the  orbit 
numbers  of  the  data  sample  are  indicated  at  the  top  of  each  graph.  Thick 
vertical  lines  are  used  to  separate  the  fit  and  predict  intervals;  thick  hori- 
zontal lines  ii’entify  the  zero-error  level.  Arrows  along  the  time  axis  denote 
the  times  of  perigee. 

Prediction  density  Models  B and  C (MODE,  MODC)  were  derived  from 
the  calorimeter  data  of  the  fit  spans.  Model  B was  based  upon  all  data  in  the 
fit  span;  Model  C was  based  on  only  the  last  orbit. 
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Fig.  A-1.  'DENSEL  (fit)/'DENSEE  (predict)  179-187 


DF:NSEL  (fit)/'DENSEL  (predict)  179-  187  (Conti 


Fig.  A-Z.  Calorimeter  (fit)/MODB  (predict)  179-187 


Fig.  A-Z.  Calorimeter  i£it)/MODR  (predict)  179-  187  (Continued) 
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Fig.  A-4.  Calorimeter  (£it)/'DENSEL,  (predict)  179-187 


i'lg.  A-4.  Calu  riniti  te  r (fit ) / 'DKNSEL  (predict)  179  187  (Continued) 
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Fig.  A-6.  Caluriiiicter  ilil)/MODB  (predict)  325-333  (Continued) 


!■ 

~ri 
1 

~f~n 

i-  ' 

T 

r 



*■ 

1 

r 

N 

rtv 

1 ; 

1 

¥ 

I ; 

■ 

■ 

{ 2 

1 1 

! ; 

_ 

• 

• 

, ; 

! 1 

1 



. 

■ 

■ 

J 

- 

i I 

t t 

1 i 

• 

• 

' i 

i i 



.—j 

’ 

( 

! i 

1 1 

I ’ 

H 

< 

_1 

- ■ • 



• - 

■ ■ 

i-  i . 

1 

_j 

. 

— 

— 

1 : 

1 

I 

■1 

■■ 

u- — 

i '. 

1 j 

■ 

BBI 

■ 

■ i ! ■ 

1 i 

■ 

■■1 

L 

, 

' 

1 T""’ 

— 

- 

mi 

1 

■n 

- i ■ 

■ 

ni 

“1 



. 

' 

I 

— 

L_ 

1 

i 

■ 

— r— 

LJ 

i! 

: 

id 

■ 

— 

H 

— 

— 

— 

— — 

tn 

— 1 — ^ — 

id 

ini 

1 

ii 

. 

. 

_J 

! ! 

T 

I 

1 

[ 

. 

L 

- 

i . - , 

1 

L . 

F 

f 

i 

1 1 ■! 

■i  ; 

I 

: 

L_ 

1 ! 

j 1 

i 

L- 

! 1 

.1 . 

1 

M 

1 ! 

! i 

w 

i 

1 , 

. ■ r 

1—  L 

r i 

■ ■ 

1 

. 

. 

' i 

.Li 

- T 

! 

3 

; 1 

1 

1 

• 

z' 

1 I 

' - 

. 

; ; ■ 

L 

1 i 

/> 

r 

•I  . ■ 

f.  -J 

r- 

( 

! 1 

■ 

r 

u 

L-L  J__ 

j ; 

J*' 

v 

T ir> 
-JO 

■ 

! i 

1 

i 

L . 

yi 

L.. 

. 

I 

«j 

i- 

i-  ! i- 

• ! i. 

I 

Cl 

• 

L_ 

LL_ 

lJ 

r 

L_ 

L 

lJ_ 

*1* 

f 4 

4 

•.  i 

_ 

aloriirK'tcr  (iit)/MC)lX,  (predict)  iZ5-5i3  ((.dnUinued  ) 


KiXS 


iSBbBmSSBBB  BSBBB^' 

^h[hbbbbb^bh 


■■■■ggl 
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-12.  Calorimeter  (fit)/'DENSEL  (predict)  50Z-510 
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Fig.  A-14.  Calorimeter  (£it)/\IODB  (predict)  7 61-769 
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Fig.  A- 15.  Calorimeter  (fit)/MODC  (predict)  761-769 
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Fig.  A-l6.  Calorimeter  (fit)/'DENSEL  (predict)  761-769  (Continued) 


Fig.  A- 17.  'DENSEL  (fit  and  predict)  792-800 
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Fig.  A-18.  Calorimeter  (fit)/MODB  (predict)  792-800 


Fig.  A-I9.  Calorimeter  (fit)/MODC  (predict)  792-800 
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Fig. 


APPENDIX  B.  Altitude  Normalization  Procedure 


The  exponential  altitude  dependence  of  the  density  generally  dominates 
over  time-  and  horizontal  space -dependent  effects,  and  thus  it  is  necessary 
to  remove  altitude  effects  in  order  to  study  atmospheric  changes.  A plot  of 
the  orbit-integrated  density  versus  altitude  of  perigee  is  shown  in  Fig.  B-1. 

A curve  of  the  form 

Y ( z ) = exp  ! c^  + c^  ( z - ) + c^  ( z - Zq 

was  fit  to  the  data  and  is  shown  in  the  figure.  In  the  expression  above  the 
c/s  are  parameters  of  the  fit,  z is  altitude,  and  z^  is  the  reference 
altitude,  chosen  to  be  70  nm. 

For  studies  of  time  and  horizontal  space  variations  in  the  drag  field, 
the  orbit-integrated  densities  { z^  ) were  reduced  to  the  reference  altitude 

Zq  using 


°o  < ‘‘o  > ' °p  ' "p 


Y ( ^0  ) / V ( .p  ) 


where  z^  was  the  perigee  altitude. 
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Fig.  B-1.  Orbit  Integrated  Density  Vs.  Altitude  of  Perigee 
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THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  V'er- 
satility  and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory;  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric optic8icKemicaTreactrons'Tn~poiluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry:  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment:  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere:  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 

THE  AEROSPACE  CORPORATION 
El  Segundo,  California 
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